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Abstract 

The study of the phenomenon of quantum weak turbulence is 
extended by determining the quasiparticle spectrum associated with 
such a system using a Green's function approach. The quasiparticle 
spectrum calculated establishes the dissipative regime and the inertial 
regime, hence a Kolmogorov type of picture. 

PACS number(s): 47.27.Ak, 67.40.Db, 67.40.Vs, 03.40.Gc 

Keywords: Kolmogorov- Zakharov spectra, Turbulence, Superflu- 
idity 



1 Introduction 

As emphasized by Newell and Zakharov the phenomena of turbulence is not 
confined simply to fluid motion. Weak turbulence in the general framework is 
the study of dispersive waves (not necessarily in fluids) with weak non-linear 
interaction. Such a system can be represented by a Hamiltonian correspond- 
ing to the wave equation with a small conservative non linearity. The study 
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of weak turbulence is relevant to many physical processes like capillary waves 
on the surface of water. 

It may be noted that the non-linear Schrodinger equation(NLS) i^!t + 
A\I/ + a| \I/ | S \I/ = 0, a candidate for weak turbulence, is widely used in 
continuum mechanics, plasma physics and optics. The NLS in the form 
i^t + A^ + a| ^ = 0, a — ±1 describes a number of physical processes 
like waves on fluid surface etc. A consequence of weak turbulence is to allow a 
Kolmogorov like spectra (Kolmogorov-Zakharov spectra) which is exact [1,2]. 
In this paper the study of weak turbulence in the quantum region described 
in Ref[3] is continued. With the help of a Green's function method the 
quasiparticle modes which describe excitations in the system are determined. 
The quasiparticle spectrum establishes clearly the dissipative regime and 
hence the range of validity of the approach. In the present study the starting 
Hamiltonian, in the language of second quantization, takes the form [3] 



The momentum conserving d dimensional 5 term has been absorbed in the T 
matrix. The above Hamiltonian may be identified as the microscopic Hamil- 
tonian for a Bose liquid. The a) , a are the usual creation and annihilation 
operators satisfying the the commutation relation [<2fc,a^] = S k i- The mode 
number operator is given by n k = a\a k - 



We interpret the state in which the occupation numbers are large as a fully 
developed turbulent system. The reason for this is that the system, as we 
will show, has a region ( inertial region) where there is no dissipation followed 
by a dissipative region at short distances. In order to study turbulence we 
replace the a k ,a\ k by c-numbers to the first approximation. This corresponds 
to assuming that the system is in a state in which there is large occupation 
number of different momenta. Such a system could model a driven system 
not in equilibrium relevant to turbulence. It may be emphasized here that 
the liquid is not a condensate which exists only at k — 0. Such a condensate 
is postulated for equilibrium calculations. We divide the Hamiltonian into 

H = H Q + i?/[3]. Since we are interested in the expectation values of 
fields between states with large mode numbers it is sensible to include as 
much of these in the unperturbed Hamiltonian H Q as the calculation permits. 
We do not need to assume that the coupling to all the quartic terms in 
the Hamiltonian is small; the diagonal part can be arbitrarily large in our 
approach. The unperturbed part (_£/<,) thus contains the quadratic part and 
the diagonal part of the T matrix i.e terms of the form T k m- After some 
simplification we can write 



Tk-^k^k^k^^ fei G^/ C2 0/ C3 (2fc 4 



(1) 




2 Turbulence in a Bose Liquid 



H Q = Y ^krik + 2^ T kt iri k ni 



(2) 
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Hi = J2 T k 1 k 2 k 3 k 4 aj k 1 ^k 2 a k . i a k4 (3) 
where we have introduced the notation T k = T kkkk ,T k i = Turn and 



rpr = f T klk2k:ik4 : ki ^ k 3 ork A ,^ 

k 1 k 2 k 3 k 4 \ q . otherwise 



The next step is to calculate (-^) , the quantum kinetic equation and look 
for the solution of the equation (^r) = 0( stationary solution). In order to 
establish Kolmogorov type solutions for the above equation the following scal- 
ing relations have been assumed : e(Xk) = X a e(k) ,T f (Xki, Xk 2 , Xk 3 , Xk^) = 
X^T'(k±, k 2 , k 3: k^ and N(e) = e z . e k is the energy corresponding to the mode 
k and N(e) is the occupation number of mode having energy e. The solutions 
corresponding to the stationary state are given by [3] 

NM(e) = C 1 e= i ^ ( 5 ) 

N^\e) = C 2 e^ (6) 

where A = M±M _ 4. The solutions mentioned above do not correspond to 
the BE distribution, hence they are the non-equilibrium solutions. We have 
ignored the two equilibrium solutions to the quantum kinetic equation which 
also exist. Having shown that such non-equilibrium scaling solutions exist 
we would like to analyse next the behaviour of the unperturbed Hamiltonian 
for the Kolmogorov type solutions as stated above. Note the unperturbed 
Hamiltonian in this approach contains the " diagonal" part of the interaction 
as described in equation (2). 



3 Energy Spectrum of the Unperturbed Hamil- 
tonian (H ) 

We are dealing with a Bose fluid, not in equilibrium, in which the occupation 
numbers of the modes are very high,i.e N k ^> 1. Hence we now proceed to 
replace the operator a k by V^fc + b k and a\ by y/Nk + 6^where N k is a 
c-number and the operators b k) tf k represent fluctuations around the ground 
state which have occupation N k for the A;th mode. This ground state will 
be represented by the ket | 0). The above substitution results in the eqn (3) 
being rewritten as 

H = Const + (£ A k b\b k + ]T B k b\b^ k + ]T B k b k b. k ) (7) 

k k k 

where 

A k = 4]T T kl Nt + u k + AT k N k (8) 
i 

B k = AT k _ k N k (9) 
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The approximate Hamiltonian,equation (7) constructed in this way can be 
called the Bogolyubov Hamiltonian. The energy spectrum for the above 
Hamiltonian may be found out by using the method of Green's functions 
as briefly described below. We define the Green's function G^\ = (0 | 
T[b{(t)b k (0)] | 0) and G^ k = (0 | T[b. k (t)b k (0)] | 0) where b\(t), &_ fc (t) are 
in Heisenberg representation. The time ordered expectation values are taken 
between the states in which there is macroscopic occupation number. This 
state is taken to represent the turbulent configuration and will be referred to 
as the ' ground state' of the system. The time derivatives of G^\ and G^ 2 \ 
yield, 

1 * k = - S (t) + 2iB k G { ~ 2 \ + iA k G^\ (10) 



dt 
dGW k 



-2iB k G<?\ - iA k G^ k (11) 



dt 

The fourier transform of the above equations with respect to the variable t 
and the evaluation of G^\uj) , the fourier transform of G^\ leads to 

G[ 1 \u) = - - ;" + 2 A * (12) 
k V ; u 2 - (A k 2 - AB k 2 ) v ; 

The poles of are given by cu = ±(A k 2 — 4i? fc 2 ) 1 ^ 2 . Hence the energy 
spectrum of the quasiparticles is given by E = E + J2k t(k)c\c k where 



e(k) = y/{A k -2B k ){A k + 2B k ) (13) 

Here E Q represents the ground state energy and c k ,c^ k are the annhilation 
and the creation operators for the quasi-particles. The equation (14) may 
be solved self-consistently using the expressions (9) and (10) and the scaling 
form of the following functions: u\ k ~ X 2 uJ k , T\ k ~ A^T^, N\ k ~ X x N k and 
e(Xk) ~ X a e(k). Let T\ k N\ k ~ X J T k N k where 7 = (3 + x. If we consider the 
term A k + 2B k , under recsaling of the momenta by A, we have, A\ k + 2B\ k = 
4A 7+1 ETmNi + AX~<T k N k + X 2 u k . In the infrared limit (k -> 0) and A -> 0, 
it is seen that the term containing uj k term is rendered ' irrelevant' for 7 < 2 
and for 7 > 2, u k is the only relevant term. Hence we have for asymptotically 
small values of k 

: 7>2 (M) 

Similar expressions can be written down for the large k spectrum. Let us 
bring in the concept of superfluidity here. Consider a body of mass M moving 
throught he fluid and in the proccess excites a quasiparticle of wave number 
k and energy e(k). Conservation of momentum yields, 

Mvi = Mv f + hk (15) 

For dissipation to occur the following relation amongst the energies should 
be satisfied, 

Mv 2 Mv f 2 

—^>—^ + e(k) (16) 
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The above two equations lead to the relation Vi > ^ where v c = ^ is 
called the critical velocity. If the velocity v of a particle moving through 
the "fluid" of quasiparticles is less than v c , where v c is the critical velocity , 
then the particle cannot loose energythe fluid is termed to be a superfluid 
i.e the external particle moving through this fluid experiences no drag. This 
is the " turbulent regime" for the system. In this regime there is no viscosity 
and hence corresponds to the large " Reynolds number" situation. The order 
of magnitude of v c is given by v c ~ ^ ~ k a ~ l . If a > 1 , it cannot be 
termed a superfluid as the critical velocity is zero. For a < 1 it is a super 
fluid as the critical velocity is very large for small k values. Let us consider 
an eddy of size k which is associated with quasiparticle of wavenumber k. 
As is seen from the expression of v c the eddy critical velocity is k dependent. 
The velocity scale of each eddy is set by a different v c . The fluid analogy is 
further strengthened by the fact that dissipation occurs when v > v c i.e for 
increasing k. We will show that a solution where the energy flow is towards 
the large k region exists. This is the turbulent flow picture. 

Note this solution makes it clear that we are in a non-equilibrium situa- 
tion. For an equilibrium configuration the system should evolve to it's lowest 
energy state(i.e small k region). 

Considering the total Hamiltonian H of equation (2), the resulting multiscale 
in the system would manifest in intermittency. It is now evident from equa- 
tion (15) that a < 1 corresponds only to the 7 < 2 case from a self-consistent 
picture. We now equate a = 7 to find out the a corresponding to the solution 
given by equations (5) and (6) . The solutions are, 

a = 3d -0 + 3 (17) 

3d- (3 = -3 (18) 

As is evident from (18), the solution (6) is unphysical. It has been explained 
later that the above fact is a consequence of the non-equilibrium nature of 
the problem. For superfluidity we have a < 1. Imposing the condition on the 
equations (16) we obtain the lower bound on (3 for the K-Z solution which 
is f3 > 11 for d=3. The result gives us the condition for superfluidity and 
also explicitly introduces a cut-off (kd) beyond which the fluid ceases to be 
superfluid(turbulent) as the quasiparticle's k dependence follows from that 
of the interaction term T^. Since we are dealing with superfluid turbulence 
here one may like to compare the energy content in the eddy k given by 
Ek = e(k)Nk to the Kolmogorov's result in the fluid turbulence , Ej, ~ k~s. 
It may be emphasized here that we are working in a region that is analogous 
to the inertial regime in the fluid turbulence and Kolmogorov's law is only 
valid in the inertial region. The infra red region in our case is defined with 
respect to the critical velocity of the system. 

For consistency of the physical picture we need to study how the energy 
flows in the model following the techniques due to Zakharov[l]. The quantum 
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kinetic equation for the mode number per unit energy range, N e , is given 
by[3], 

{ dt ' 



J de 1 de 2 de 3 [U(e 1 ,e 2 ,e 3 ,e)(e 1 e 2 e 3 e) z (e 3 z + e z - 
(e y + e 3 y - e 2 y - e^)5(e + e 3 - e 2 - e l )} (19) 

where y = 3z — p — 3,N e ~ e z and p = M±M _ 4 Tj nc [ er the rescaling of the 
variables e,e 1 ,e 2 and e 3 by A we have U(Xei, Xe 2 , Ae 3 , Ae) = A p f/(ei, e 2 , e 3 , e) 
where p = 3d+2 P — 4. Denoting the above integral by K(t, z, y) we have the 
following relationi^(e, z, y) = A _!/_1 i i C(Ae, z, y). Putting A = K Hence we 
have 

K(e,z,y) = e- v - 1 K(l,z,y) (20) 

We now define the current corresponding to the conservation equation for 
the mode number per unit energy range(iV e ) as j. The conservation equation 
is, 

<^> = ! = ^*(i.„) (2D 

^3 = -— K(l,z,y) (22) 

y 

j can be identified with the energy curent per unit energy range as can be 
seen from simple dimensional analysis. Considering the stationary solutions 
i.e. K-Z solutions we have for the case y — 0, 



dK(l,v,y) 

J = — 



o (using L'Hospital's rule) (23) 



dy 

Considering the other stationary solution i.e y =1 we have, 

j = (24) 

Since we are considering turbulence, a non-equilibrium phenomenon, only 
the solution corresponding to j 7^ is the true solution while j = is the 
unphysical one. Thus the flux analysis rules out the second solution. The 
fact that a cannot be self-consistently determined from the second solution 
is manifest in the unphysical nature of the solution as discussed above. Note 
that the current corresponding to the non-zero solution of equation (21) 
represents an energy flux from low k end (wavnenumber at which energy is 
injected into the system) to the high k end where dissipation takes place in a 
manner very similar to one encounters in ordinary fluid turbulence governed 
by the Navier-Stokes equation (See Fig 1). 

The results confirm that only one of the scaling solutions found corre- 
spond to a non- equilibrium situation in which the energy flux is non-zero. 
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Superfluid Viscous 



k = k — > oo 

Direction of energy flow 



Figl 



In a realistic situation there will be a driving source term for injecting en- 
ergy into the system. Such systems have been extensively analysed by Za- 
kharov[l]. It is clear that the arguments continue to hold for the general 
quantum system we have considered. The analysis of the quasiparticle spec- 
trum associated with quantum turbulence thus clearly establishes the inertial 
range, it places restrictions on parameters, demonstrates that scaling solu- 
tions which are not in equilibrium are possible and for the case a < 1 i.e where 
v c is k dependent suggests that a novel form of intermittency can occur. The 
analysis thus completes the approach outlined in Ref[3]. We conclude with 
the following remarks. In the general approach of Zakharov to weak wave 
turbulence it is shown [1] that non-equilibrium scaling solutions are possible 
and that there is flux of either momentum or energy present in the system. 
Our earlier analysis of the quantum version of Zakharov approach [3] estab- 
lished that these features are also present in the quantum system described. 
In this paper we have gone further. It has been shown that underlying such 
a non-equilibrium solution there is an inertial range with no dissipation and 
that a self consistently determined short distance dissipative scale can be 
found. Thus for the quantum system considered the analysis shows that the 
original intuition of the Kolmogorov approach to scaling can be justified. 
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